


The HEGRA Experiment 
Status and Recent Results 

F.A. Aharonian^ and G. Hcinzclmann ^ for the HEGRA collaboration 

'^Max-Planck-Institut ffir Kcrnphysik, P.O.BOX 103980, D-69029 Heidelberg, Germany 

''Umversitat Hamburg, H. Inst, fiir Exp. Physik, Lumper Chaussee 149, D-22761 Hamburg, Germany 



tin 



> 

o 

(N 

O 

a^ 

:^ 

Oh' 

o 

U 



X 



A report is given on the status and recent results from the HEGRA imaging atmospheric Cherenkov telescopes 
(lACT) concerning Mrk 501 and GRS 1915+105 and on recent results from the HEGRA arrays concerning galactic 
and extragalactic 7-sources, counterparts of TeV Gamma-Ray-Bursts, cosmic ray anisotropies, energy spectrum 



and composition of primary cosmic rays. 
Cosmic Ray Symposium 1996. 



The report summarizes the presentations given at the 15 European 



1. GENERAL 

The HEGRA (High Energy Gamma Ray Astron- 
omy) experiment is located at the Observatorio 
del Roque de los Muchachos on the Canary Island 
of La Palma, (2200 m a.s.L, 28.75° N, 17.89" W) 
covering an area of about 200m • 200m. The lay- 
out is shown in Fig.^. The experiment consists of 
6 imaging air Cherenkov telescopes and 3 arrays, 
the scintillator array, the wide angle Cherenkov 
array AIROBICC ||l| and the Geiger tower ar- 
ray[|. 

One of the major aims of the HEGRA ex- 
periment is the search for galactic and extra- 
galactic sources of TeV 7-radiation and a de- 
tailed study of their characteristics with the imag- 
ing atmospheric Cherenkov technique at energies 
above 500 GeV. Hereby the stereoscopic obser- 
vation with the system of telescopes will be ex- 
ploited. The TeV 7-astronomy is continued to- 
wards higher energies (> 20 TeV) with the three 
arrays. For the arrays a further important aim is 
a measurement of the energy spectrum and ele- 
mental composition of primary cosmic rays in the 
energy region from « 500 TeV to w 10 PeV. 

In the following first the status and perfor- 
mance of the imaging atmospheric Cherenkov 
telescopes and their recent results concerning 
Mrk 501 and the microquasar GRS1915-I-105 are 
given, then the recent results from the arrays are 
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243 Scintillation Counters (1 m ) 
77 AIROBICC Detectors ( I, 0.4 m) 



6 Cherenkov Telescopes 
17 Geiger Towers (18m") 



Figure 1. Layout of the HEGRA experiment 

presented concerning galactic and extragalactic 
7 sources, TeV counterparts of GRB's, cosmic 
ray anisotropies and the energy spectrum and el- 
emental composition of the primary cosmic rays. 



2. THE HEGRA SYSTEM OF lACTs 



2.1. Telescopes 

Presently we are very close to coniplete a stereo- 
scopic system consisting of 5 imaging telescopes 
(reported by G. Hermann and M.Ulrich at this 
symposium) which during the next several years 
will serve as one of the most powerful instruments 
in the field, operating at energies above 500 GeV 

i- 

The multi mirror reflector of each telescope 
(Fig. H) has an area S'mir — 8.5 m^. The 60 cm di- 
ameter front aluminized and quartz coated spher- 
ical glass mirrors with a focal length of about 5 m 
are independently mounted on an almost spher- 
ical frame of an alt-azimuth mount. The angu- 
lar resolution of the reflector is better than 10 
arcminutes. The mount is driven by two stepping 
motors, and controlled by shaft encoders which 
provide a tracking accuracy of 0.02°. 

Each telescope is equipped with a 271-channel 
camera [Q consisting of a close-packed cluster of 
circular PMTs (EMI 9083 KAFL) faced with hol- 
low hexagonal-to-round light guides (Fig ||) . The 
diameter of a pixel is about 0.24°, which results 
in the telescope's field of view FoV« 4.6°. The 
photon to photoelectron conversion efficiency in- 
cluding all kinds of losses (due to the mirror re- 
flectivity, the quantum efficiency of PMTs, losses 
in light guides, etc.), and averaged over the spec- 
trum of the Cherenkov light at the level of obser- 
vations, is around jyph^e = 0.15. 

The entire electronics, including the trigger, 
readout, and the telescope control, is based on 
the VME standard. The digitization of the PMT 
pulses is performed with a Flash ADC system 
which allows to reconstruct the shape and time 
distribution of the pulses. The trigger decision is 
taken on two levels [^. The fast trigger requires 
a minimum number of pixels to be fired, and the 
second (slower) level restricts the topology of an 
image. Near each telescope there is a small hut 
which houses the local electronics. Each telescope 
is equipped with a local CPU which communi- 
cates via an Ethernet line with the central com- 
puter. Although each telescope is built as stand 
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Figure 2. The HEGRA telescope CT-3. 



alone detector, all telescopes can be integrated 
into the system by software command through 
the central computer (Pentium PC). 

The first nominal 'system-telescope' called (for 
chronological reasons) CT-3, is in operation since 
August 1995. The data taken in the 95/96 sea- 
son generally confirm the expected performance 
of a single telescope. The next 3 telescopes CT-4, 
CT-5, and CT-6 have been installed, tested and 
calibrated during July-December 1996, and are 
now operating together with CT-3 as a system 
in the 2-fold coincidence mode. The preliminary 
results are very encouraging and already demon- 
strate the power of the stereoscopic approach (see 
below). Currently, we use different FADC sys- 
tems in different telescopes, but in the near future 
(until March 1997) we plan to operate the above 
four telescopes with the 120 MHz FADC system. 
At approximately the same time we plan to re- 
place the present 61-channel (pixel size 0.43°, and 
FoV~ 3.9°) camera of CT-2 by a standard 271- 
pixel high resolution camera. Thus we hope that 
after the test and calibration work, we will be able 
to run in the second half of 1997 the whole system 
as a very sensitive instrument for deep study of 
VHE 7-rays in the energy region 0.5-10 TeV. 

In addition we plan to upgrade the first (pro- 
totype) HEGRA telescope CT-1. This equato- 
rial mount telescope with 5 m^ reflector was in- 
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Figure 3. The 271-channel camera. 



stalled and is taking data since 1992 [g. In 1994 
the original 37-channel camera was replaced by 
a new high resolution 127-channel camera (pixel 
size 0.24°) based on the same pixels used in the 
'system' telescopes, but designed in the CAMAC 
standard 0. In 1997 the round glass mirrors will 
be replaced by light hexagonal aluminum mirrors 
with a total reflector area ~ 12 m^. Thus the en- 
ergy threshold of the telescope will be reduced to 
500 GeV. We plan to use this telescope for indi- 
vidual observational programs, especially in im- 
mediate reaction to potential VHE 7-ray emitters 
showing flaring activity in other wavelengths. At 
the same time this telescope being located very 
near to the central telescope CT-3, can be effec- 
tively integrated, if necessary, with other lACTs 
in order to set up 2 or 3 independent stereoscopic 
lACT systems for simultaneous observations of 
different sources. 

2.2. Performance of the lACT system 

Although 7-ray showers of energy 500 GeV with 
core position within 100 m from the telescope pro- 
duce > 100 photoelectrons in the camera, suffi- 
cient for proper image analysis, the energy thresh- 
old of a single telescope is somewhat larger, and 
is determined by the suppression efficiency of the 
night sky background (N.S.B.) at the trigger level. 



Figure 4. The hardware differential detection rates 
of cosmic rays and 7-rays by CT-3. The rates for 7- 
rayes are calculated assuming a power-law spectrum 
with spectral indices: (1) - q=3; (2) - a— 2.5; (3) - 
a=2, and for an absolute fiux of 10"^^ ph/cm^s above 
1 TeV. 



The hardware trigger criterion for each telescope, 
optimized by Monte Carlo calculations, is orga- 
nized such that the signals in at least 2 neigh- 
boring pixels from the inner 169 pixels should ex- 
ceed ~ 15 photoelectrons (2W/169 > 15). The 
Monte Carlo calculations of the detection rates 
of 7-rays from a point source, presented in Fig. 
m show that the effective energy threshold of the 
telescope CT-3 is between 700 GeV and 1 TeV, 
depending on the slope of the primary photon 
spectrum. The first observations with CT-3 gen- 
erally confirm the calculated performance of a sin- 
gle telescope, namely, a detection rate of 7-rays 
from the Crab Nebula above 700 GeV of about 
1.5 X 10^"^ photons/s on top of a cosmic ray (CR) 
detection rate « 5 Hz, the significance of the sig- 
nal per 1 hour of observations being ~ 2.2(t. 

The requirement of simultaneous detection of 
showers by locally triggered > 2 telescopes signif- 
icantly reduces the effect of the N.S.B. as well as 
that of the another important background con- 
nected with local muons. This allows a lower en- 
ergy threshold of the system down to < 500 GeV. 
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Figure 6. The reconstructed coordinates of the 
coincidence events observed by CT2 and CT3 within 
0.3° of the direction to the Crab (ON source events) 
and offset of the Crab (OFF source events) for 8.5 
hour observations in both cases. 



Figure 5. ALPHA-distribution for the Crab in the 
camera of CT-3 after the application of 'standard' 
image cuts in > 2 telescopes. 



Also, the ability of the system to analyze the 
Cherenkov images in different projections where 
the images are only partially correlated, provides 
a superior rejection of hadronic showers. The pre- 
liminary results of observations of the Crab Neb- 
ula in October '96 by CT3/CT4/CT5, operating 
in the 2-fold telescope coincidence mode, gener- 
ally confirm this prediction. In Fig. o we present 
the ON-source and OFF-source ALPHA distri- 
butions for the Crab in the camera of CT-3, ob- 
tained after filtering the events through the 'stan- 
dard' image cuts of the showers detected by > 2 
telescopes. This preliminary result shows that 
the significance of detection of ~ 500 GeV 7-rays 
from the Crab is achieved at a level of ~ Scr per 
1 hour. In fact, we expect further significant im- 
provement of the signal using energy and distance 
dependent 'cuts' to the detected images. (Note 
that the stereoscopic approach allows to deter- 
mine the energy and the distance to the shower 
core with accuracies < 20% and ^ 10 m, respec- 
tively). Furthermore, the obvious advantage of 
the joint operation of > 2 telescopes is the pos- 



sibility to determine, without any a priory as- 
sumption, the arrival directions of individual 7- 
ray primaries with accuracy < 0.1°. For the case 
of point-like sources, this provides by itself, re- 
jection of the CR background based only on the 
orientational characteristics of showers by a factor 
of at least several 100. The procedure obviously 
requires good understanding of the characteristics 
of all telescopes linked in the system. Although 
the work regarding the CT3/CT4/CT5 data set is 
still in progress, our results obtained earlier using 
the combination of telescopes CT1/CT2 ^, and 
telescopes CT2/CT3 allows us an optimistic view 
concerning the practical realization of this impor- 
tant aspect of the stereoscopic approach. Indeed, 
despite the short time of simultaneous observa- 
tions from the Crab by CT-2 and CT-3 taken dur- 
ing the season 1995/1996 and a low coincidence 
rate (due to the high energy threshold of CT-2), 
the stereoscopic approach reveals a clear signal 
with 23 'excess' events against 4 (!) 'background' 
events within an angle 0.3° around the direction 
to the Crab (see Fig. g), the latter being in good 
agreement with the Monte Carlo predictions for 
the angular resolution of the CT2/CT3 system. 

Thus, our current understanding of the system 
based on the recent experimental results as well 
as on comprehensive Monte Carlo studies, allows 
us to make conclusive predictions for the perfor- 



mance of the system of the 5 HEGRA imaging 
telescopes, namely 

• energy threshold around 500 GeV; 

• angular and energy resolutions of about 0.1° 
and < 20%, respectively; 

• minimum detectable fluxes of 7-rays .Fmin ^ 
10^^-^ ph/cm^s from point sources, and Fmin ~ 
10^^^ ph/cm^s from moderately extended sources 
with an angular size ~ 1°. 



2.3. The recent results 

While the 5-IACT system is nearing the final 
stage of its commission, we continue to use the 
available telescopes in different configurations not 
only for methodological purposes, but also for ob- 
servations of a number of astronomical objects. 
The first positive signal of 7-rays was seen from 
the Crab Nebula in 1992 by CT-1 in its first ver- 
sion M . Later we have detected TeV 7-rays from 
the Crab with a significance of ~ lOcr using CT-2 
l9| . The current sensitivity of the system consist- 
ing of CT3/CT4/CT5 has achieved a level which 
allows to obtain a strong signal from the Crab 
for just a few hours of observations under almost 
background- free conditions (see Fig. 0). Thus, 
we can now use the high energy photon beam of 
this persistent source for calibration purposes and 
for an understanding of detector characteristics. 
Over the last several years CT-1 and CT-2 have 
been used for the search for TcV 7-rays from other 
objects, in particular Cygnus X-3, AE Aquary, 
etc., reported earlier by different groups as possi- 
ble episodic 7-ray emitters (see e.g. |1C[| ). Unfor- 
tunately, we did not find DC 7-ray signals from 
any one of these sources. However, the obser- 
vations of Markarian 421 in 1995 resulted in the 
independent detection of statistically significant 
positive signals by CT-1 and CT-2 ^^, thus con- 
firming the discovery by the Whipple collabora- 
tion of TeV emission from this BL Lac object |12[ . 

In 1996 we allocated an essential part of the 
time of the available telescopes for observations 
of relatively nearby AGNs, in particular Mrk 421 
(CTl, CT2, CT3), Mrk 501, PKS 2209, MS 
01166 (CTl), W Comae (CTl, CT2, and CT3, in 
the framework of the multiwavelength campaign). 



Table 1 

Summary of Mrk 501 and Crab observations with 
CT-1. 



object 


Mrk 501 


Crab 


ON-source obs. time [h] 


146.8 


22.9 


ON events after all cuts 


1325 


423 


expected background 


974 


237 


number of excess events 


351 


186 


excess rate [h~^] 


2.4 


8.1 


stat. significance (Li & Ma) 


5.2 


7.6 



and 1ES2344+514 (CT3/CT4/CT5). Also, CT2 
and CT3 were used for observations of the su- 
perluminal galactic source GRS 1915-1-105, and 
CT3/CT4/CT5 for observations of the shell type 
supernova remnant 7 Cygni. The analysis of most 
of these data is still in progress. Here we briefly 
discuss the results obtained from Mrk 501 and 
GRS 1915+105. 

Mrk 501 (reported by S. Bradbury) 

This source is very similar to Mrk 421, and rep- 
resents the same sub-population of AGNs called 
BL Lac objects - highly variable AGNs with- 
out strong emission lines, but showing a strong 
nonthermal component of radiation from radio 
to X-ray wavelengths. Recent EGRET observa- 
tions showed that these objects are effective 7-ray 
emitters at energies > 100 MeV |13|. Although 
Mrk 501 is not in the list of 7-ray emitting BL 
Lac's detected by EGRET, this source has been 
reported by the Whipple group to radiate in 7- 
rays above 300 GeV at an average flux level of 
approximately 0.1 Crab [Q. The low flux of 
VHE 7-radiation from this source makes its detec- 
tion rather difficult for our CT-1 and CT-2 tele- 
scopes. However, at the beginning of this year 
Mrk 501 started to show a fiaring activity at VHE 
7-rays (Weekes, private communication). This 
motivated us to allocate the CT-1 for long-term 
observations of Mrk 501. In order to obtain max- 
imum exposure time the source was observed in 
tracking mode. Almost 150 hours of good quality 
data was obtained at zenith angles < 25° during 
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Figure 7. The ALPHA-distribution of the excess 
events for Mrk 501. 



the period March- August 1996. The OFF-source 
observations required for background determina- 
tion were made when Mrk 501 was not observable. 

The ALPHA distribution for Mrk 501 shown in 
Fig. ^ indicates an excess of 351 events against 
expected 974 ± 57 background events (see Ta- 
ble El) which can be interpreted as detection of 
7-rays with a rate R^ ~ 2.4 h~^ and statistical 
significance 5.2ct |i3|. Using the same procedure 
we find an excess rate of 7-like events from the 
direction of the Crab, R^ ~ 8.1 h^^. This im- 
plies that during the observation period the flux 
from Mrk 501 was « 3 times lower than the Crab 
flux, i.e. approximately 2.5 x lO^^^ph/cm^s. For 
comparison, the Whipple observations of Mrk 501 
in 1995 showed lower flux relative to the Crab 
(« 10%) at energies > 300 GeV. This can be ex- 
plained by different energy spectra of the Crab 
and Mrk 501 and/or by a strong variability of 
Mrk 501 which can be significant on timescales of 
1 day [Q. The daily excess 7-ray event rates of 
our data are shown in Fig. ra. Above our thresh- 
old of 1.5 TeV we have seen no significant "burst 
behaviour" . Although the analysis of the distri- 
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Figure 8. 
Mrk 501. 



Daily excess 7-ray event rates from 



bution of total light content in the camera does 
not show a significant deviation of the spectral 
shape of Mrk 501 from the Crab spectrum, due 
to the uncertainty in the estimated power-law in- 
dex 2.6 ± 0.5 we cannot make strong statements 
about modifications to the power law spectrum 
such as cut offs. 

GRS 1915+105 (reported by F.A. Aharonian) 

The recent discovery of the galactic super- 
luminal source GRS 1915-fl05 |Q, the scaled 
down analog of the AGN jets which are recog- 
nized as high energy 7-ray emitters, motivated 
the HEGRA collaboration to put this "micro- 
quasar" into the "Targets of Opportunity" list. 

The first observations of the source in Spring 
1995 with the telescopes CTl (13 h) and CT2 
(25 h) in ON/OFF mode showed a small evidence 
for a signal with statistical significance 2.8ct and 
2.6(7, respectively. However, the observations by 
the same telescopes after Aug 15 with doubled 
exposition time did not reveal any indication of 
a signal. Since this might be a result of a highly 
sporadic behavior of GRS 1915-f 105, we decided 
to continue observations in 1996. 

In 1996 GRS 1915-^105 has been observed dur- 
ing the period May-July 1996 by CT2 and CT3. 

CT2 with an energy threshold « 1 TeV had 
been successfully used earlier for observations 
of the Crab Igll and Mrk 421 ill. The main 




Table 2 
Summary of GRS 1915+105 observations with 
CT-2. 

May June July Total 
ON-source obs. 

time [h] 9.89 37.12 23.83 70.84 

ON events after 

all cuts 427 1527 944 2898 

expected 

background 365 1373 875 2612 

number of 

excess events 62 154 69 285 

excess rate [li^i] 6.2 4.2 2.9 4.02 

stat. significance 
(Li & Ma) 2.68 3.46 1.95 4.66 



Figure 9. The ALPHA-distribution of the excess 
events for GRS 1915+105. 



characteristics of this telescope are studied by 
Monte Carlo simulations which reproduce quite 
well the distributions of different image parame- 
ters. The rejection power of the telescope based 
on both shape and orientation parameters allows 
to suppress the CR background by a factor of 
100, the 7-ray acceptance being at the level of 
> 50%. The sensitivity of the telescope based 
on the observations of the Crab Nebula is es- 
timated as « 1.25cr/h detection rate of a sig- 
nal from a point source at the flux level J-y{> 
ITeV) ~ 10""'^^ ph/cm^s. Remarkably, this pa- 
rameter has not been changed noticeably during 
the last 2 years which demonstrates the high sta- 
bility of operation of the telescope. 

Since GRS 1915+105 is a highly variable ob- 
ject, we observed the source in the tracking mode. 
In addition, each month certain ON/OFF pairs 
were taken in order to monitor the shape of the 
background "a-distribution" . The statistics of 
observations at zenith angles 9 < 30° and the 
results are summarized in Table ^. The ALPHA- 
distribution of ON-OFF source events shown in 
Fig. g indicate the existence of a possible signal 



with significance 4.7cr. The average rate of excess 
events is ~ 4h~^ which corresponds to w 25% 
of the 7-ray detection from the Crab. Note that 
most of the "excess" events are due to the June 
period of observations which can be explained by 
~ 50% of statistics contained in this period. At 
the same time, perhaps one cannot exclude an 
intrinsic variability of the signal itself. For ex- 
ample the statistical significance of the signal in 
June 19 alone is of about 3cr, with a rate exceed- 
ing the average rate by a factor of 3. This is seen 
in Fig. ^ where we present the daily rates of 
the "excess" events. Another feature of the time 
profile of the rate is its rather interesting coin- 
cidence with 3 prominent few-day-duration fiares 
of GRS 1915+105 detected in radio and X-rays 
(see Fig. |o|) . This possible correlation in time of 
the rate of "excess" events may serve as an addi- 
tional argument that strengthens the belief that 
the TeV 7-ray detection from this source is a real 
effect rather than a result of statistical fluctua- 
tion of the background. If so, assuming that the 
radiation originates in the jet of GRS 1915+105 
approaching us, the intrinsic average TeV lumi- 
nosity of the source is estimated as 10'^^ erg/s. 

However, we would like to emphasize that be- 
fore going too far in the astrophysical interpreta- 
tion of this statistically marginal signal, we obvi- 
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Figure 10. Daily rates of the excess events from 
GRS 1915+105. For comparison, the periods of large 
outbursts in radio and X-rays also are shown. The 
light hatched zones (1,2,3) with borders indicated by 
dashed vertical lines corresponds to the periods of 
X-ray flares reported by Greiner J. et al., (Ap.J. Let- 
ters (in press); lAU Circ.6435), and the dark hatched 
zones (1*,2*,3*) correspond to the radio flares de- 
tected by the Ryle and Nancay telescopes (f AU Circ. 
6411; 6432). A strong hard X-ray outburst in May 
24-27 (MJD 50228-50231) has been reported also by 
the BATSE team (lAU Circ.6411). Note that at the 
time intervals without 7-ray data points correspond 
to the moonlight periods when the source could not 
be monitored. 



ously need independent confirmations of this re- 
sult by other instruments and groups. As was 
mentioned before, the source was observed si- 
multaneously also by the CT-3 telescope. At 
the Perpignan Symposium, a positive signal was 
also shown from this telescope. Further analysis 
of these data however has revealed evidence for 
hardware problems during this observation period 
in some of the camera pixels, resulting in large 
spurious signals, which may fake 7-ray topologies. 
Work is in progress to correct these problems in 
software and to reprocess the data. At this point, 
no statements concerning GRS 1915-1-105 signals 
in CT-3 data can be made. 



3. RESULTS FROM THE ARRAYS 

3.1. Angular resolution of AIROBICC 

In an analysis using Monte Carlo and experimen- 
tal data it has been shown that the angular reso- 
lution of the AIROBICC array above a threshold 
of 18 TeV for all events inside the array border 
is (763% = 0.29 °. It decreases to 0.1 ° at higher 
energies. The corresponding angular resolution 
of the scintillator array above 20 TeV is fTg3% = 
0.9° and reaches at higher energies 0.4 ° ||l^. The 
numbers given refer to hadron-initiated showers 
of vertical incidence. For photon-initiated show- 
ers an angular resolution of cr63% — 0.22 ° is ex- 
pected according to Monte Carlo calculations for 
the AIROBICC array. 

3.2. Separation of gamma- and hadron- 
initiated showers with pearl 

A method to achieve a gamma-hadron separation 
has been presented by J. Prahl. The separation of 
gamma- and hadron-initiated showers is, besides 
a good angular resolution, crucial for a search of 
7-ray sources. 

The method called pearl uses first the well 
know fact that electromagnetic showers decline 
after the shower maximum faster than hadronic 
showers. Therefore the ratio No/Cl of parti- 
cles measured at ground Ng and the number of 
Cherenkov photons Cl is a well known separa- 
tion quantity pq ]. 

The other newly introduced quantity in the 
method pearl is related to the larger transverse 
extent of hadronic showers as compared to elec- 
tromagnetic showers due to the hadronic trans- 
verse momentum arising in hadronic interactions. 
The quantity used for the transverse spread is de- 
fined by rp ■■=J2i\\'^i-'^c\\ rir I Y.i'^i whereby 
Tc is the reconstructed core position, f^ the po- 
sition of the i* scintillator counter and n^ the 
number of particles in the i scintillator counter. 
The separation quantity rp is related to the age 
parameter from a NKG fit. 

Instead of using quantities obtained from a fit 
to the lateral distributions the method pearl uses 
the corresponding quantities derived from count- 
ing the amplitudes in the different counters i.e. 
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Figure 11. The separation quantities rp vs. 

N(e)/N(L) (s. text) for photon- and hadron-initiated 
showers (Monte Carlo). The straight line (pearl cut) 
indicated accepts 48 % photon showers and rejects 
90 % hadron showers (E > 18 TeV). 



the sum of the pulseheights in the scintillators 
N(e),the sum of the Cherenkovlight measured in 
the AIROBICC array N(L), and the quantity rp 
sensitive to the lateral extent defined above. The 
advantage of these quantities is that showers at 
the trigger threshold can be analysed without re- 
quiring a fit to converge. (Studies have shown 
that for showers for which the lateral fits con- 
verge the results are the same in the energy range 
of interest). 

In Fig. ni^ the two separation quantities are 
plotted against each other for photon and hadron 
initiated showers. A clear distinction can be seen 
and a straight line (pearl cut) can be used to sep- 
arate photon- and hadron initiated shower in the 
2-dim logarithmic plot. The Q factor defined as 
obtained for this method varies 



from Q=1.5 at threshold of 18 TeV to Q=2.5 at 
50 TeV. 

3.3. Galactic TeV 7-ray sources 

With the gamma-hadron separation method 
pearl described above, based on the scintilla- 
tor and AIROBICC arrays, a search for galac- 
tic sources has been performed using the 94/95 
data. Amongst 13 selected sources the most sig- 
nificant excess of 3.3 a after gamma-hadron sep- 
aration is observed for Her X-1. No excess has 
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Figure 12. Integral 7-ray flux for the Crab 

nebula from the HEGRA Arrays (scintillator and 
AIROBICC) together with results from the HEGRA 
Cherenkov telescopes and other experiments. 



been observed for the crab nebula in the 93/95 
data whereas in the 1992/93 data a 3.6 a effect 
has been reported (24*/i ICRC Rom 1995, un- 
published). The results for the Crab nebula are 
summarized in Fig. |l^. 

3.4. Extragalactic TeV 7-ray sources 

Up to now two extragalactic TeV 7 sources Mrk 
421 and Mrk 501 have been observed with the 
technique of imaging air Cherenkov telescopes 
(lACT) up to < 10 TeV. 

Whereas with the lACT only one source at a 
time can be observed, the arrays allow a simulta- 
neous observation of several sources in the field of 
view of HEGRA, however for the prize of a higher 
energy threshold. This may severely limit the 
chance of an observation due to the absorption 
of TeV photons on the not well known infrared- 
optical background fields (7-ray horizon). 

In order to perform a search a list of 11 source 
candidates with redshifts < 0.1 has been set up, 
belonging to the class of blazars visible in the field 
of view of HEGRA [M. The sources have been 
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looked at in the stacking mode and individually. 
The search has been performed for datasets 
with and without gamma-hadron separation ca- 
pability. The gamma-hadron separation with the 
GEIGER towers using a neural net analysis (Q- 
factors 2 - 2.5) ||2^] requires a threshold of 50 
TeV and an angular bin of 1 degree (using the 
direction of the scintillator array) therefore the 
results from other data sets are for comparison 
also evaluated for 50 TeV and 1 degree. Stack- 
ing all source candidates the highest significances 
reported by N. Magnussen are observed in the 
scintillator data from 6.89 -6.92 without gamma- 
hadron separation (3.5 a) and from 11.94-3.95 
with gamma-hadron separation using the Geiger 
towers (4.0 a). The most prominent single source 
candidate is 0116+319 which contributes with 
4.4 a in the scintillator data. No excess has been 
observed for a sample of blazars with higher red- 
shifts and for a sample of 49 galaxies character- 
ized by steep radio spectra. 



As reported by J. Prahl the scintillator data 
from 1994/95 show an excess of 1.2 a and the 
AIROBICC data from 1994/95 with 7-hadron 
separation 2.4 a. Exploiting the much better an- 
gular resolution and lower energy thresholds pos- 
sible for the AIROBICC data, the significances 
do not increase as expected but decrease. In or- 
der to compare the results of the different data 
sets the excesses with the highest significance are 
converted tentatively into a flux and the other 
results in upper hmits (Fig. p^ ). 

A definite conclusion cannot yet be drawn and 
more systematic checks are under way especially 
concerning the absolute pointing which is crucial 
in order to exploit the good angular resolution of 
the AIROBICC array. Also further AIROBICC 
data from 1992/93 are investigated. If a signal 
could be established this would imply, that the 
infrared background radiation has to be much less 
than presently assumed and the sources according 
to present day models would probably be proton 
accelerators. 
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Figure 13. Comparison of the results from different 
analyses and data samples (s. text) on a flux scale for 
a generic blazar (from a superposition of 11 nearby 
blazars) 



3.5. TeV counterparts of gamma-ray bursts 

Gamma-Ray-Bursts (GRB's) are known for more 
than 25 years yet their physical origin and dis- 
tance scale are still unknown. Observations of 
TeV counterparts would impose limits on the un- 
known distance scale because of the absorption at 
the infrared intergalactic background fields which 
limits the mean free path of TeV Photons to 
distances< 100 Mpc. 

To search for TeV emission associated with 
gamma-ray bursts registered by the BATSE in- 
strument all HEGRA detector components have 
been used. Especially searches with the large field 
of view of the air shower arrays have been per- 
formed for TeV emission which may occur simul- 
taneous, preceding or following a burst. A de- 
layed emission may be expected if GRB's account 
for the acceleration of protons up to the highest 
energies which in turn generate a secondary in- 
tergalactic electromagnetic cascade due to inter- 
actions with the microwave background fields. 

No significant signal has been observed for the 
different searches as reported by B. Funk and H. 
Krawczynski. The results from a search for co- 
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Figure 14. For 31 BATSE bursts (1994-1996) the 
HEGRA upper limits on the coincident integral flux 
(90% Confidence Level) as a function of the threshold 
energy for inducing photons are shown. The thresh- 
old energy is a function of the zenith angle of the 
location of a burst. 
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Figure 15. ^-y/^CR upper limit (90% CL) from 
the galactic disc (S=scintillator, A=AIROBICC, 7= 
with pearl 7-hadron separation). Numbers are Li-Ma 
significances. 



incident emission obtained for 31 BATSE bursts 
are shown in Fig. ^ To guide the eye an extrap- 
olated GRB spectrum |2^ is displayed. 

In future a lower energy threshold of about 2 
TeV will be possible due to higher data taking 
rate and an angle sensitive trigger. 

3.6. Cosmic ray anisotropies 

A search for anisotropies has been reported by 
D. Schmele using the data taken during 1994/95 
with the scintillator array (60 Mio events ), with 
the AIROBICC array (50 Mio) and in the data 
with the AIROBICC array using the 7-hadron 
separation method pearl. 

In order to search for anisotropies, all data (for 
each data set separately) have been used to deter- 
mine a two dimensional shape of the background 
distribution in local coordinates cj) and 6. This 
shape of the background distribution is corrected 
for the air pressure by a multiplicative term de- 
pending on 6 and air pressure. Using this (pres- 
sure corrected) background shape in local coordi- 
nates a background map in equatorial coordinates 
is generated according to the actual event rates 
recorded. In order to obtain the actual event rates 
in this procedure, the experimental event rates 
have first been corrected for the (variable) dead- 
time and air pressure. In this way the actual up 



time and performance of the array is taken into 
account. 

The deviations of the data from the background 
(for bin sizes corresponding to the angular reso- 
lution) follow a normal distribution N(0,1) rather 
well. For binsizes from 1° to 50° an upper limit 
of ^cb/^bg < 10~^ for the smallest bin and 
< 10~^ for the largest bin is obtained. These val- 
ues are essentially determined by the statistics of 
^ 50Mio events per event sample. 

In order to look for an emission from the galac- 
tic disk, the maps are transformed into galactic 
coordinates and projected onto galactic latitudes 
and compared. A disk of width ||6|| < 10° has 
been assumed (1 = 30° - 220°). The results for 
the upper limits for a flux from the galactic disk 
for the three data sets are shown in Fig.^. Also 
shown are the limits of CASA MIA J2|], the ex- 
pected theoretical flux of 3x 10~^ for diffuse emis- 
sion from interaction of cosmic rays with inter- 
stellar matter ||6|| < 5° H and a naive EGRET 
extrapolation with a constant slope of 7 — 2. The 
HEGRA limits of ^■y/<^CR > 4 x lO""* can be in- 
terpreted as limits for an additional flux due to 
'unresolved' galactic sources in the energy range 
above 20 TeV. 

Finally a search for a contribution from the 
region around the direction of the EHE events 
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as seen by Fly's Eye and Yakutsk has been per- 
formed. A contribution in the TeV energy range 
may be expected via the intergalactic cascade 
{p + 1S.K -^ X A- jTev)- On the scale from 5° 
to 20° radius around the Fly's Eye event position 
(RA=85.2, DEC=48.0) ||, no excess has been 
found. An upper limit for the contribution of an 
EHE source ^source/^BG ^ 10"^ above 20 TeV 
can be given. This limit includes a possible con- 
tribution from the galactic disk, which is near to 
the estimated source position. 

All searches are carried out with the three data 
samples described above and all three agree well 
within their different sensitivities and statistics. 
In the near future the size of the data sample will 
will be increased substantially. 

3.7. Primary energy and composition 

The measurement of the energy and composition 
of primary cosmic rays in the energy range of the 
knee is only accessible to the indirect methods of 
air showers and is still controversial and a mat- 
ter of debate. Here a new method is presented 
to reconstruct the energy and masses of primary 
particles with energies above 300 TeV from the 
electromagnetic component of an air shower as 
measured by a scintillator array and an array of 
Cherenkov counters. The main steps of the recon- 
struction procedure as presented by H. Krawczyn- 
ski, can be summarized as follows: 

(1) The distance to the shower maximum D„iax 
is reconstructed from the lateral distribution of 
the Cherenkov light as measured at detection 
level. It can be shown by Monte Carlo sim- 
ulation that Dniax can be determined indepen- 
dent of the primary mass (and almost indepen- 
dent of energy). This is due to the fact that 
the shower development behind the shower max- 
imum, which dominates the determination of the 
distance to the maximum, is independent of the 
primary mass of nuclei. 

(2) Electromagnetic energy Eem- From the 
number of electrons Ne as obtained from a NKG 
fit and the distance to the showermaximum D„iax 
the number of electrons in the maximum N^^max 
can be inferred. This again is independent of the 
the particle mass. The number of electrons in 
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Figure 16. The upper part shows the distribution of 
the reconstructed energies. On the lower part the dis- 
tribution of the reconstructed nucleon number of all 
events with reconstructed energies between 500TeV 
and IPeV is shown. The abundances of three groups 
of primaries have been determined with a three com- 
ponent fit of the distributions predicted from Monte 
Carlo data for each group (dotted lines) to the mea- 
sured distribution (solid line). 



the maximum is in a good approximation pro- 
portional to the electromagnetic energy Eem- 

(3) Energy per nucleon E/A. From the distance 
to the showermaximum the penetration depth is 
obtained. Since the penetration depth depends 
(logarithmically) on E/A, the energy per nucleon 
E/A can be (coarsely) determined. 

(4) Primary energy E. The fraction of the 
total energy which is transformed into electro- 
magnetic energy increases in a known way (i.e. 
from Monte Carlo simulations) with the en- 
ergy/nuleon: Ee„i=f(E/A). With E/A from step 
(3) and Ee„i from step (2) the total energy E can 
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be determined. 

(5) Primary nucleon number A. With E and 
E/A known from step (3) and (4) also the nucleon 
number A is coarsely known. More details of the 
method can be found in Pq| . 

In order to first test the new method, it is ap- 
plied to the energy range from 500 TeV to 1 PeV 
where the composition is known or can be ex- 
trapolated from direct measurements taken with 
satellites and balloone-borne experiment. The 
preliminary results shown in Fig. nfl are in good 
agreement with the results from direct measure- 
ments. For the differential energy spectrum a 
spectral index of 2.60±0.10±0.1 is obtained in the 
energy range from 500 TeV to 2 PeV. The abun- 
dance of light nuclei to all nuclei is obtained to 
be 0.62±0.05±0.10 in the energy range from 500 
TeV to 1 PeV. The distributions in Fig.[6| reflect 
the fluctuations in the penetration depth for the 
different primaries. It should be stressed that for 
the mass determination the mean values as well as 
width (i.e. fluctuations) of the curves contribute. 

Results with a slightly different method have 
been reported at 24*'' ICRC Rom ^ . In the near 
future the new method is applied to the energy 
range above 1 PeV (knee). 
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